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ABSTRACT 


T'.'.e  cbieccive  or  this  paper  is  to  develop  accepcanee  procedures  rased  on  a 
sequential  probability  ratio  test  when  the  system  failures  follow  Weibull  proba¬ 
bility  distribution.  The  acceptance  procedures  in  both  Milicarv  Handbook  103  an 
Mil-Sta-731C  are  based  on  exponential  distribution  of  failures.  Exponential  dis 
tribucion  implies  a  constanc  hazard  rate.  The  Veibull  distribution  is  much  mere 
comprehensive  and,  with  suitable  choice  of  parameters,  can  describe  constant, 
increasing  or  decreasing  hazard  rates. 

The  specif ications  are  usually  written  in  terms  of  the  scale  parameter  : ,  c 
shape  parameter  d  being  determined  by  the  nature  of  the  system  under  test.  Mien 
N’  units  are  being  tested  and  n  units  have  failed  at  time  t,  the  acceptance  is 
based  on  the  monitoring  of  a  statistic  :r  given  as 

:  .  •  ~n  ( log  a .  -  lov  q  ')  f  — ■  7  c'  +  (N-r.K  ‘  * 

"  °  _i=l  ^ 

wnera  o  ,  and  ,i  ,  are  the  desirable  and  unacceocab ie  values  of  the  scale  oarameter 

1  i 

rhe  indicated  decision  is  to 

accept  if  i  a  constant,  log  Z 
reject  if  :r  i  another  constant,  log  Z, 

and  ti  continue  testing  if  remains  within  the  two  limits. 


SEQUENTIAL  TESTS  UNDER  WEI BULL  DISTRIBUTION 


I .  Objective 

The  objective  of  this  paper  is  the  development  of  acceptance  procedures 
based  on  a  sequential  probability  ratio  test  when  the  system  failures  follow 
a  Weibull  probability  distribution.  The  acceptance  procedures  in  both  Military- 
Handbook  H  108  and  standard  Mil-Std-781C  are  based  on  an  exponential  distribu¬ 
tions  of  failures. 

There  are  two  major  reasons  for  attempting  to  base  acceptance  procedures 
on  a  Weibull  distribution  rather  than  on  an  exponential  distribution. 

1.  Weibull  distribution  is  more  comprehensive  than  exponential. 
Exponential  distribution  would  only  apply  to  those  systems 
where  the  hazard  rate  remains  constant  over  time.  Weibull 
distribution  would,  with  suitable  choice  of  parameters,  be 
applicable  to  all  systems  where  the  hazard  race  is  constant, 
increasing,  or  decreasing. 

2.  The  modei  describing  the  growth  of  the  reliability  of  a  new 
system  under  development  can  be  shown  to  be  a  special  case  of 
the  Weibull  model.  The  growth  model,  first  studied  by  Duane 
(1),  describes  a  sequence  of  changing  system  conf igurat ions  and 
is  different  from  the  static  reliability  models  for  a  fixed 
configuration.  However  the  model  itself  follows  a  probability 
distribution  which  can  be  approximated  by  Weibull  distribution 
with  shape  parameter  d  =  -1'2. 

’.'he  principal  reason  for  a  sequential  probability  ratio  test  is  that 
such  test  procedures  involve  a  minimum  amount  of  testing. 


II.  Probability  Ratio  Tests  for  Mult i -parameter  Functions 


A.  Svmbols 


Let  f(x,  6.,  0^,  '5j_)  be  probability  density  function  where 


a  discrete  valued  variable 


t^i’  =  parameters 


=  a  c-componertt  vector  0^,  6 

■=  a  k-component  vector  $  ,  .... 

=  parameter  space  for  parameters  0^,  . 


• 0,. 


iij,  =>  parameter  space  for  parameters  4h 

a,b  *  subspaces  of  and 

I  s. 

n 

L  -  II  f(x.,T,K) 


L*(a,b)  =  max  L|Tca,  Keb 

\(a,b)  =  probability  ratio 

L* (a ,b) 

L*(nT,s„) 

1  fv 


B.  General  Probability  Ratio  Test 

The  parameter  set  T  is  the  set  under  stud  .  Set  K  includes  all  other 
parameters.  Hence  all  hypotheses  are  stated  in  terms  of  ".  There  will  be 
3  cases  of  the  tests  of  hypotheses  according  to  the  assumption  about  K. 

Under  assumptions  of  f ( • )  satisfying  general  regularity  conditions  and 
large  n,  in  a  test  of  the  hypothesis 

V  1  ’  T° 

the  quantity  -2  Log  is  approximately  distributed  as  chi-square  with 

t  degrees  at  freedom  when  is  true. 

This  approximate  distribution  can  he  used  to  devise  a  test  of  the 
hypothes i s 


1 


against  a  composite  hypothesis 
H^:  T  /  T° 

The  power  of  this  test  is  not  easy  to  determine  as  the  distribution  of 

-2  log  )  is  non-central  chi-square. 

K. 

C.  Special  Cases  of  Probability  Ratio  Tests 

Three  special  cases  of  probability  ratio  tests  can  be  defined  according 
to  the  assumptions  regarding  K. 

Case  1  K  is  assumed  to  be  known. 

Case  2  An  estimate  K  is  used  in  place  of  K. 

Case  3  K  is  unspecified. 

Case  1:  This  case  involves  t lie  testing  of  a  simple  hypothesis 


agains^  a  simple  hypothesis 
H  :  T  =  T1 

witii  an  added  condition  that  K  is  known  to  be  equai  to  a  constant  K.  In  this 
case  the  parameter  space  under  both  Hq  and  would  be  single  points  and 

L* ( • ,  ■ )  =  L 

For  this  case  a  test,  witii  specified  values  x  and  2  of  Type  I  and  I'vpe  il 
errors,  can  be  devised  following  Wald  (.-i)  . 

If  is  true  and  i  false,  then 

L(T°,n)  -  l-i 


-3- 


and 


UT1,^)  <  d 


Therefore  the  test  statistic 


satisfies 


s 


1 


L(T°.K)  1-a 

UT1,?)  " 


i 


If  Hq  is  false  and  is  true,  then 
L(T°,K)  <  a 

and 

KT^K)  ^  i_e 

Therefore  tiie  test  statistic  s ^  satisfies 

_<  '  i 

’L  -  L-i 

Inequalities  (1.1)  and  (.J.i)  define  the.  acceptance  and  rejection  regions 

1  1—1 

resoectivelv .  When  s,  lies  between  — and  — —  there  is  no  decision.  in 
i  1-d  a 

this  case,  the  sample  site  can  be  sequentially  increased  until  an  accept  or 
reject  decision  is  readied.  Such  a  test  would  be  called  a  Sequential 
Probability  Katie  E'est,  SPRT. 

Sequential  Probability  Ratio  Test  was  first  developed  by  A.  W'alu  y  ■-  i  . 
The  optimality  o>‘  SPRT  was  formally  proved  by  A.  Wald  uul  .  Wolfowic:  '  "d  . 
The  problem  of  •  >d  t  ima  1  i  t  v  of  SPKT  for  Markov  processes  was  studied  by  ah.osii 


(2).  A.  Sirjaev  (3)  considers  the  more  general  optimal  stopping  problem  for 


Markov  processes.  Sirjaev  (3)  also  gives  a  discrete  time  formulation  of 
SPRT. 

Case  2:  This  case  also  would  be  like  Case  1  but  with  an  assumed  value 
K  being  used  instead  of  the  true  value  K..  In  this  case  the  computed  values 
of  L  would  be  in  error  and  these  errors  can  be  defined  as  error  factors. 


e0  " 


_  L(T°,K) 


L(T°,K) 


and 


=  l(tl,k) 

1  L( X1 ,K) 


Then  the  test  statistic  s^  would  be  given  b\ 


and  it  would  be  possible  to  define  the  acceptance  and  rejection  regions  of 
tiie  tvpe 


accept 

if 

S3  -  q 

reject 

if 

s3  :  q 

and  increase  sample  size  if  q^  <  <  q^. 

An  approximating  distribution  of  the  statistic  s,  would  cave  to  be 
developed  to  obtain  the  Type  I  and  Type  II  errors  for  t.iis  test. 


-o- 
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Case  3:  This  case  involves  the  casting  of  a  simple  hypothesis 


V  T  -  T 


against  a  simple  hypothesis 


H1 :  T  =  T 


Consider  the  test  statistic  s0  given  by 


MT°,\)  l*it°,.ok)  L*(«T,nK) 


S2  “  1  •  L*(QT,i\) 


x(T',aK) 


T  ’  K  L*(T*,»y 


If  H  is  true  then  L*(T°,Q  )  >  1-a  and 
U  K.  — 


L*(T1,l'J  )  <  B. 


Consequently 


1-a 

S2  >  ~ 


Similarly  when  is  true 


2  -  1-3 


Again,  as  in  the  cases  (1)  and  (2)  these  two  inequalities  define  die 
acceptance  and  rejection  regions.  When  s^  lies  between  and  there 


is  no  decision  and  the  sample  size  can  be  sequentially  increased  until  an 
recept  or  reject  decision  is  reached. 


III.  Reliability  Tests  Using  Welbull  Distribution 

A.  Weibull  distribution  is  a  3  parameter  distribution  with  parameters 
q  =  scale  parameter 
6  =  shape  parameter 
and  y  =  location  parameter 

For  reliability  studies  the  location  parameter  is  known.  In  "time-to-failure” 
tests  y  is  always  zero.  In  the  "stress-to-f ailure"  tests  y  is  known  from 
engineering  considerations.  Then  the  Weibull  distribution  can  be  transformed 
to  a  two  parameter  distribution  by  counting  the  time  (scress)  as  t-y  instead 
of  t. 

Further  the  reliability  specifications  are  written  in  terms  of  MTTF . 

For  a  specified  value  of  8  these  can  be  seen  as  specifications  for  q.  The 
Table  6-1  in  the  appendix  gives  the  values  of  a  corresponding  to  the  values 
of  MTTF  for  specified  9. 

A  good  estimate  of  9  is  almost  always  available  from  historic  data  about 
similar  projects.  In  the  proposed  acceptance  procedure  the  historic  estimate 
of  8  can  be  used  in  place  of  the  known  9  and  this  assumption  can  be  validated 
by  testing  an  associate  hypothesis  H^:  S  =  9  at  frequent  intervals  during 
the  testing. 

Although  it  is  assumed  here  that  a  historic  estimate  of  9  is  available, 
there  can  be  circumstances  where  such  estimates  of  9  are  not  available.  For 
example  a  totally  new  system  may  not  have  suitable  data  available.  In  these 
cases  the  analysis  would  follow  "case  3"  in  the  preceding  section.  It  would 

be  difficult  to  obtain  any  operational  procedures,  however,  as  the  computation 

0  0  - 
of  L*(T  would  not  be  as  straightforward  as  the  computation  of  L(T  ,K). 

K. 

This  case  is  not  addressed  to  in  the  following  operational  procedure. 


-7- 


B.  Symbols 


qo: 

An  acceptable  value  of  q 

V 

An  unacceptable  value  of 

<1 

6: 

Historic  shape  parameter 

estimate 

t.  : 
l 

Time  to  the  iC^  failure 

t: 

Time 

N: 

Number  of  units  being  tested 

n: 

Number  of  units  failed 

C.  Test  of  Che  Hypothesis 

The  test  of  Che  hypothesis 


V  q  =  q0 


against 


Hr  q  ■ 


under  the  assumption  of  9  being  known  would  come  under  Case  1.  Here  ?  =  {q 
and  K  =■  {9}.  For  the  stipulated  Type  I  and  Type  II  error  probabilities  of 
and  a  respectively,  the  test  itself  will  take  the  form 

„  .  ,  1— a 

accept  if  s.  >  — — 
i  “  b 

reject  if  sL  <  y— 


,  .  „ .  .  --  a  1-a 

ana  continue  testing  if  -j—g  <  <  — — 


L(T°,K)  -  (  n  q0t.“Exp(-q0t.S+1/9+l)]-Exp(-(N-n)q0te4‘1/e+l) 
i*l  L 


-8- 


and 


1—0  0+1  0+1 
L(T  ,K).  =  [  IT  q^t  Exp (-q  ^t  ^  /'c+1)  ]  ■  Exp  ( IN— n ) q  ^ t  '/0+1.) 

i=l 


Taking  logarithms  and  simplifying 


,  ,  L(T° ,K) 

log  s  =  log  - r2^- 

L(TA,K) 


=  n( log  aQ  -  log  o^) 


qrqo  r  n  0+1  ,  0+1, 

~0+l~  , “  ti  +  (N-U)t  J 

i=l 


Let  z  =  -i— i  and  z  =  — ii-_  .  Than  the  acceptance  condition  becomes 
Ip  Z  1-p 


log  s,  ^  log 


and  the  rejection  condition  becomes 


log  a  <  log  z. 


Further  algebraic  simplification  leads  to  the  conditions  based  on  another 
statistic 


s*(t)  =  Z  t,9+1  +  (N-n) t 9+1 
l-l  1 


so  that 


<1,  -  <3  . 

L  0 

log  s  ^  =  -n(  log  q  ^  -  log  q())  +  +  y  s*,t) 


leading  to 


accept  if  o  >  log 
reject  if  o  <  log  z7 

and  continue  sampling  if  o  remains  within  the  two  limits. 


D.  Computational  Procedure 

It  can  be  noticed  that  the  right  hand  sides  of  the  two  conditions  are 
constant.  These  can  be  plotted  as  two  horizontal  lines  and  :  can  he  p  . ..  itou 
as  a  function  of  time. 

At  any  interval  when  there  is  no  failure  o  would  be  decreasing  and 
there  would  be  sharp  upward  jumps  whenever  there  is  a  failure.  The  test 
chart  would  appear  as  given  in  Figure  (i). 

It  can  be  further  observed  that  when  there  is  no  failure,  the  term 

1  ”  ^  0 

-n(log  -  log  q^)  remains  constant  and  the  term  ^  —  1 —  s*(t)  in  an 
increasing  function  of  time.  The  slope  of  tiiis  function  at  time  t  can  be 
calculated  as 


0  +  1 


s*(t) j 


e  +  l 


d_ 

dt 


s*(t) 


and 


d_ 

dt 


s*(t) 


9+1 


t .  + 


i 


(N-n) 


td+i] 


=  (N-n) (9+1) t ' 


-10- 


Figure  1 


Sequential  Acceptance  Chart 


^  1  ~  a 

slope  =  — ; - : —  (N-n)(0+l)t" 

b  +  1 


=  (c; x  -  qQ)  (N-n)t9 


At  every  failure  epoc  the  first  term  would  decrease  by  the 
log  qL  -  log  q0). 

lo  facilitate  the  computations  the  quantities  t  have  been 


quantity 


tabulated 


or  selected  values  of  9  ill  fable  (b-2)  in  Lhe  appendix. 


IV .  Operational  Procedure 

The  operational  procedure  for  sequential  tests  under  Weibull  distribution 
can  be  summarized  as  follows. 

1)  Determine  the  desired  and  unacceptable  values  of  MTTF. 

2)  Obtain  an  estimate  of  8.  Preferably  choose  an  acceptable  8  value 
out  of  the  following 

-1,  -0.5,  0,  0.5,  1,  2,  3. 

Almost  any  practical  failure  density  can  be  closely  approximated 
by  one  of  these  values. 

3)  Obtain  q^  and  q^  from  Table  (6.1). 

4)  Determine  N  the  number  of  units  to  be  tested. 

5)  Maintain  and  update  the  chart,  similar  to  Figure  (1). 

Mote:  When  a  large  number  of  units  are  to  be  tested  it  would  be  desirable 
to  revalidate  the  estimate  of  0  at  suitable  intervals. 

Appendix 

Table  6.1  Values  of  q  for  stated  MTTF 

0 

Table  6.2  Values  of  t 


VI. 
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Values  of  Parameter  q  for  Selected  MTTF  Value 


m 

S  -  -0. 

50 

9  =  0.00 

e  =  l.oo 

6  =  2.00 

e  =  3. 

00 

1 

0.  707 

0.100 

1 

0.157 

1 

0.213 

1 

0.269 

1 

2 

0.500 

0.500 

0.392 

0.267 

0.168 

3 

0.408 

0.333 

0.174 

0.791 

-1 

0.333 

-1 

4 

0.353 

0.250 

0.981 

-1 

0.333 

-1 

0.105 

-1 

5 

0.316 

0.200 

0.628 

-1 

0.170 

-1 

0.431 

<■> 

6 

0.288 

0.166 

0.436 

-1 

0.988 

-2 

0.208 

_2 

7 

0.267 

0.142 

0.320 

-1 

0.622 

_2 

0.112 

-2 

8 

0.250 

0.125 

0.245 

-1 

0.417 

_2 

0.659 

-3 

9 

0.235 

0.111 

0.193 

-1 

0.293 

-2 

0.411 

-3 

10 

0.223 

0.100 

0.157 

-1 

0.213 

-2 

0.269 

-3 

15 

0.182 

0.666 

-1 

0.698 

-2 

0.632 

-3 

0.533 

-4 

20 

0.158 

0.500 

-1 

0.392 

^  L. 

0.267 

-3 

0.168 

-4 

25 

0.141 

0.400 

-1 

0.251 

-2 

0. 136 

-3 

0.691 

-5 

30 

0.129 

0.333 

-1 

0.174 

-2 

0.791 

-4 

0.333 

-5 

35 

0.119 

0.285 

-1 

0.128 

-2 

0.498 

-4 

0.179 

_ 

40 

0.111 

0.250 

-1 

0.981 

-3 

0.333 

-4 

0.105 

-5 

45 

0.105 

0.222 

-1 

0.775 

-3 

0.234 

-4 

0.658 

-6 

50 

0.100 

0.200 

-1 

0.628 

-3 

0.170 

-4 

0.431 

-6 

55 

0.953 

-1 

0.181 

-1 

0.519 

-3 

0.128 

-4 

0.295 

-6 

60 

0.912 

-1 

0.166 

-1 

0.436 

-3 

0.988 

-5 

0.208 

-6 

65 

0.877 

-1 

0.153 

-1 

0.371 

-3 

0.777 

-5 

0.151 

-6 

70 

0.845 

-1 

0.142 

-1 

0.320 

-3 

0.622 

-5 

0.112 

-6 

75 

0.816 

-1 

0.133 

-1 

0.279 

-3 

0.506 

-5 

0.853 

-7 

80 

0.790 

-1 

0.125 

-1 

0.245 

-3 

0.417 

-5 

0.659 

-7 

85 

0.766 

-1 

0.117 

-1 

0.217 

-3 

0.347 

-5 

0.517 

-7 

90 

0.745 

-1 

0.111 

-1 

0.193 

-3 

0.293 

-5 

0.411 

-7 

95 

0.725 

-1 

0.105 

-1 

0.174 

-3 

0.249 

-5 

0.331 

-7 

100 

0.707 

-1 

0.100 

-1 

0.157 

-3 

0.213 

-5 

0.269 

-7 

200 

0.500 

-1 

0.500 

_2 

0.392 

-4 

0.267 

-6 

0.168 

-8 

300 

0.408 

-1 

0.333 

_2 

0.174 

-4 

0.791 

-7 

0.333 

_o 

400 

0.353 

-1 

0.250 

-2 

0.981 

-5 

0.333 

-7 

0.105 

-9 

500 

0.316 

-1 

0.200 

0.628 

-5 

0.170 

-7 

0.431 

-10 

600 

0.288 

-1 

0.166 

_2 

0.436 

-5 

0.988 

-8 

0.208 

-10 

700 

0.267 

-1 

0.142 

-2 

0.320 

-5 

0.622 

-8 

0.112 

-10 

800 

0.250 

-1 

0.125 

-2 

0.245 

-5 

0.417 

-8 

0.659 

-11 

900 

0.235 

-1 

0.111 

-2 

0.193 

-5 

0.293 

-8 

0.411 

- 1 1 

1000 

0.223 

-1 

0.100 

-2 

0.157 

-5 

0.213 

-8 

0.269 

-11 

2000 

0.158 

-1 

0.500 

-3 

0.392 

-6 

0.267 

-9 

0.168 

-12 

3000 

0.129 

-1 

0.333 

-3 

0.174 

-6 

0.791 

-10 

0.333 

-13 

4000 

0.111 

-1 

0.250 

-3 

0.981 

-7 

0.333 

-10 

0.105 

-13 

5000 

0.100 

-1 

0.200 

-3 

0.628 

-7 

0.170 

-10 

0.431 

-14 

6000 

0.912 

-2 

0.166 

-3 

0.436 

-7 

0.988 

-11 

0.208 

-14 

7000 

0.845 

-2 

0.142 

-3 

0.320 

-7 

0.622 

-11 

0.112 

-14 

8000 

0.  790 

-2 

0.125 

-3 

0.245 

-7 

0.417 

-11 

0 . 65q 

-15 

9000 

0.745 

-2 

0.111 

-3 

0.193 

-7 

0.293 

-11 

0.411 

-15 

10000 

0.707 

-2 

0.100 

-3 

0.157 

-7 

0.213 

-11 

0.269 

-15 

Mote : 

In  this 

tab  Le 

the  values 

of  q 

are  listed  in 

the  form 

x*lo"" 

For 

txamp le 

j  0. 125 { 2 

means  0.125*102 

or 

12.5. 

Table  6.2 


Values  of 


a 

t 


t 

9  =  -0.50 

6  =  0.00 

9  =  1.00 

9  =  2.00 

6  =  3.00 

1 

0.100 

1 

0.100 

1 

0.100 

1 

0.100 

1 

0.100 

l 

1 

0.707 

0.100 

1 

0.200 

1 

0.400 

1 

0.800 

l 

3 

0.577 

0.100 

1 

0.300 

1 

0.900 

1 

0.270 

•> 

4 

0.500 

0.100 

1 

0.400 

1 

0.160 

2 

0.640 

2 

5 

0 . 447 

0.100 

1 

0.500 

1 

0.250 

O 

0.125 

3 

6 

0.408 

0.100 

1 

0.600 

1 

0.360 

2 

0.216 

3 

7 

0.377 

0.100 

1 

0.700 

1 

0.490 

0 

0.343 

3 

3 

0.353 

0.100 

1 

0.800 

1 

0.640 

2 

0.512 

3 

9 

0.333 

0.100 

1 

0.900 

1 

0.810 

2 

0.729 

3 

10 

0.316 

0.100 

1 

0.100 

2 

0.100 

3 

0. 100 

4 

15 

0.258 

0.100 

1 

0.150 

2 

0.225 

3 

0.337 

4 

20 

0.223 

0.100 

1 

0.200 

2 

0.400 

3 

0.800 

4 

25 

0.200 

0.100 

1 

0.250 

2 

0.625 

3 

0. 156 

5 

30 

0.182 

0.100 

1 

0.300 

2 

0.900 

3 

0.270 

3 

35 

0.169 

0.100 

1 

0.350 

2 

0.122 

>4 

0.42S 

5 

40 

0.158 

0.100 

1 

0,400 

2 

0.160 

4 

0.6*0 

5 

45 

0.149 

0.100 

1 

0.450 

2 

0.202 

4 

0.911 

5 

50 

0.141 

0.100 

1 

0.500 

0.250 

4 

0.125 

0 

55 

0.134 

0.100 

1 

0.550 

2 

0.302 

4 

0.2  66 

b 

60 

0.129 

0.100 

1 

0.600 

> 

0.360 

4 

0.216 

b 

65 

0.124 

0.100 

1 

0.650 

2 

0.422 

4 

0.274 

b 

70 

0.119 

0.100 

1 

0.700 

2 

0.490 

4 

0.343 

b 

75 

0.115 

0.100 

1 

0.750 

7 

0.562 

4 

0.421 

b 

80 

0.111 

0.100 

1 

0.800 

■7 

0.640 

4 

0.512 

6 

35 

0.108 

0.100 

1 

0.850 

•7 

0.722 

4 

0.614 

6 

90 

0.105 

0.100 

1 

0.900 

•7 

0.310 

4 

0.  '29 

o 

95 

0.102 

0.100 

1 

0.950 

7 

0.902 

4 

0.857 

b 

100 

0.100 

0.100 

1 

0. 100 

3 

0.100 

5 

0 . 100 

; 

200 

0.707 

-1 

0.100 

1 

0.200 

3 

0 . 400 

5 

0 . 800 

7 

300 

0.577 

-1 

0.100 

1 

0.300 

3 

0.^00 

5 

0 . 2  70 

S 

400 

0.500 

-1 

0.100 

1 

0.400 

3 

0.  lbC 

b 

0 . 640 

s 

500 

0.447 

-1 

0.100 

1 

0.500 

3 

0.250 

b 

'7.125 

■~4 

oOO 

0.408 

-1 

0 . 100 

1 

0 . 600 

3 

0.360 

n 

0 . 2  lo 

700 

0.377 

-1 

0.100 

1 

0.700 

3 

0.490 

b 

0.  3*3 

*4 

800 

0.353 

-1 

0.100 

1 

0.800 

3 

0 . 640 

0 

0.512 

9 

900 

0.333 

-1 

0.100 

1 

0.900 

3 

0.810 

b 

0.729 

Q 

1000 

0.316 

-1 

0.100 

1 

0.100 

4 

0.100 

7 

0 . 100 

10 

2000 

0.223 

-1 

0.100 

1 

0 . 200 

400 

7 

0 . 800 

10 

3000 

0. 182 

-1 

0.100 

1 

0.300 

4 

0 . 11 00 

; 

0 . 2  70 

’  i 

4000 

0.153 

-1 

0.100 

1 

0 .  *00 

-4 

4  .  i6() 

3 

1 1 .  fi40 

1  1 

3000 

0.141 

-1 

0.100 

1 

0.5<I0 

*4 

0 . 250 

s 

i).  125 

l: 

6000 

0.  i29 

-1 

0.100 

1 

0 . 600 

^4 

0. 360 

s 

■> .  2  i  <> 

.  j 

7000 

0.119 

-1 

0.100 

1 

0.700 

9.  *90 

0 

O.  3*3 

1 3 

3000 

0.111 

1 

—  4. 

0. 100 

1 

0.800 

>4 

0 .  '-.*() 

8 

0  .  i  1. 

i  j 

9000 

0.105 

-1 

0 . 100 

1 

0.900 

•4 

0,  .  ■{  [i.! 

6 

1 1 .  7  2  i 

i: 

10000 

0.100 

-1 

0.100 

1 

0.100 

3 

0.100 

) 

'i.  1 1 ; • ' 
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